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Abstract 

We propose a new estimate of the process B — > (p, oj) + 7 by using the light-cone 
QCD sum rule. The aim is to choose the chiral quark current operators in order to 
eliminate the uncertainty arising from the chiral-even operators in the correlator. In the 
case of neglecting the corrections from p- meson mass, the sum rules for form factors are 
obtained, which are valid to the twist-three accuracy. The resulting branching ratios 
are consistent with the experimental data. 
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The rare decays of B mesons have been studied for years for a precise understanding of 
the physics of CP- Violation related to the structure of the CKM mixing matrix elements 
and testing of standard model to search for new physics. At present, only upper bounds on 
the exclusive modes B — > (p, uo) + 7 are available from CLEOII, namely Br(B° — > p 0/ ~f) < 
3.9 x KT 5 , Br(B° -> ^7) < 1.3 x 10~ 5 , and Br(B~ -> p~-y) < 1.1 xlO -5 in Ref.fl]. The rare 
decay modes B — > (p, u) + 7 are very suitable for study with the light-cone QCD sum rule 
method, which is developed in Ref.[2,3] recently, and have been investigated in Refs.[4,5]. The 
points of the light-cone QCD sum rule are combining the traditional QCD sum rule method 
with the hadronic wavefunction descriptions of the hard exclusive process. However, instead 
of the small distance x ~ in the traditional QCD sum rule method, the operator product 
expansion (OPE) in the light-cone sum rule is carried out in terms of nonlocal operators near 
the light cone x 2 ~ 0; and instead of the vacuum condensate, the nonperturbative dynamics 
is parametrized as so-called light-cone wavefunctions classified by their twist. Using the light- 
cone QCD sum rule, the authors in Ref. [4] found the form factor F B ^ P (0) = 0.285 ± 15% 
with the b-quark pole mass m b = 4.8GeV and a threshold s = 33.5 ± 0.5GeV 2 . However, 
g\ and g^_, the nonleading twist distribution wavefunctions of p meson used in Ref. [4,5], are 
not known well with a drastic effect on the resulting precision. In addition, there are the 
uncertainties in the b-quark pole mass and condensate parameters, which is an important 
input parameter in the light-cone sum rule for the form factor F B ^ P (0). Therefore, a further 
precision investigation of this rare process is valuable and necessary. Ref. [6] has made a new 
estimate of the process B — > K*j by choosing the chiral current operators in the correlator. 
Similarly, through the choices of the chiral current operators in the correlation function, the 
uncertain nonperturbative quantities can be eliminated effectively from the respective sum 
rules for the form factor F B ^ P (0) and decay constant f B , and thus a more precise prediction 
for the rare process can be made. 

Following the procedure in Ref. [6], we calculate the rare decay B — > pj which is induced 
by the effective Hamiltonian 
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where £ 4 = V t bV t * d , are Wilson coefficients and the apparent forms of operators Oi(fi) 

are discussed in Refs.[7,8]. 

In the process B — > p7, the electromagnetic penguin operator 7 = -^^m^da^ ^-^-bF^ 
plays the most important role and it leads to the effective Hamiltonian responsible to b — > 

d + r- 

H eff = Cm^da^il + 75 )g"6, (2) 

where e and q are the emitted photon polarization vector and momentum, respectively; the 
constant C depends on the CKM matrix elements V t bV t * d . 
Then the amplitude for B(p + q) — > pQo, 77) + 7(5, e) becomes: 

A(S^p7) = Crn^e" < p(p, 77)1^(1 + l5 )q ,/ b\B(p + q) > . (3) 

With cr M1 , 75 = |e M zy a/ 30" a/3 , we can parametrize the hadronic matrix element (at q 2 = 0) in 
Eq.(3) as 

< p{p, V )\da^(l + l b )q v b\B{p + q) >= (-2ie^7/"gV + 2p • q^ ~ 2? • t^F^O). (4) 

Thus the decay width T(B — > /ry) can be easily written as 

r(B - fry) = ^G 2 F \V th f \V t *f |F^(0)f C 7 ( Xt , m h f{ml + m 2 d ) { ^^ . (5) 

Now we are in a position to calculate the form factor F B ^ P (0) in the light-cone QCD sum 
rule framework. In order to eliminate the uncertain nonperturbative quantities g v L and g", 
we adopt the following correlator of the chiral current operators 

Fvfaq) = * J d A xe^ x < p(p,r])\Td(x)a^(l + ^)q v b(x)b(0)i(l + l5 )u(0)\0 >, (6) 

which is different from Ref. [4,5]. F^{p,q) can further be written as two separate terms: 

F^(p,q) = ijd^xe^ <p(p,77)|T(d(x)(7^(l + 75)g , '6(a;)6(0)m(0))|0> 

+ < p(p, r])\T(d(x)a^(l + l5 )q»b(x)b(0)i l5 u(0))\0 > . (7) 

Both pseudoscalar resonance states B^ and scalar resonance states B^ in the hadronic ex- 
pression give contributions to the correlation function: 

< p\Tido^{l + i b )<fb)\B >< B\bi l5 u\0 > 



F^(p,q) = 



m 2 B -(p + q) 2 
2 



< p\T{da, v {\ + ^b)\B h p{s) >< B h p{s) \U{\ + l5 )u\0 > 
+ 2- —2 ~uTT7a2 W 



h m ln - (P + ?) 2 

p(s) 



which can be parametrized as the Lorentz-covariant form: 

F^(p, g) = (-2ie^ Q/3 *f gV + 2p • - 2g • VP ,)F((p + g) 2 ). (9) 

The invariant amplitude F((p + q) 2 ) obeys a general dispersion relation in the momentum 
squared {p + g) 2 : 

F((p + q) 2 ) = / — ds + subtraction terms. (10) 

■/m B s - (p + q) 2 

Using the relation < Blbi^ulO >= —fs, and isolating the contribution of the lowest B 
meson from the spectral density, we have 



m. 



P(s) = S(s - ml)-*f B F B ^(0) + p h (s), (11) 
m b 

where p h (s) denotes the spectral density of higher resonances and continuum states B^, 
and is approximated as by invoking so-called quark-hadron duality ansatz 

p h (s)=p QCD (s)Q(s-s ) (12) 

with the perturbative spectral density Pqcd{ s ) and the threshold parameter s Q . Putting 
everything together, we obtain the hadronic expression for F((p + q) 2 ): 

On the other hand, the form factor F^((p + q) 2 ) is calculable in the region of large spacelike 
momenta (p + q) 2 <C by the OPE near the light cone x 2 = 0. Because of the large virtuality 
of the b quark, it is safe to calculate the form factor in the light-cone expansion. Here we 
confine ourselves to the leading contribution in the OPE, which comes from the contraction 
of the b-quark operators to the free b-quark propagator: 

< 0|TK*)S(0)|0 > = / *Eje-*-i£±^ (14) 



Iq lQn 2 a 2 2a' 



It has been shown in Ref.[4] that the gluon radiative corrections make less than 5% contri- 
bution and can be ignored in our calculation. Substituting Eq.(14) into Eq.(6), we obtain 

F,(p,q) = 2m b jf -^-^ j dxe^~ m ^ < v ) \d(x)a^(l + 75 )«(0)|0 > . (15) 

With the correlator we choose, only the chiral current matrix element < p(p, rf)\d{x)a IJlv q u (\ + 
7 5 )m(0)|0 > remains, while the chiral-even operator matrix elements < p|J(x)7 M w(0)|0 > and 

< p|g?(x)7 m 75-u(0)|0 > disappear. This will be very useful to improve the sum rule for the 
form factor F B ^ P (0). If we neglect the corrections from p-meson mass, the matrix element 

< p(p, rj)\d(x)a^ u q u (l + ^ 5 )u(0)\B > can be defined as Ref.[9] 

< p{p,ri)\d{x)<T IJa ,^{l+'y 5 )u(0)\0 >= i[{q-r])p ll -{p-q)r]^ie lxva pe v q a pf 3 \ f due iup - x f±<p±(u), 

(16) 

where the variable u is the fraction of p-meson momentum carried by the d quark and the 
value of decay constant f± can be found in Ref.[10], f± = (160 ± 10) MeV; the leading twist 
wavefunction (p±(u) should be taken at the scale characterizing virtuality of the b quark, 
p 2 ~ m% - ml ~ 5GeV 2 . 
At present, the F((p + q) 2 ) in light-cone QCD sum rule reads 

F((p + q) 2 ) = m b / — — / dxe^-<^- a / due iu ^ ± (u, p 2 )/ ± . (17) 
Jo Lb7r z a J Jo 

Completing the integral over x and a, and replacing the variable u by s = ^ + (1 - u)m 2 p , 
we have the form factor 

F((p + g) 2 ) = m b r ^f\ 2 2 " 2 f±ds. (18) 
Jm\ s — (p + q) 2 u 2 m 2 p + m b 

Making the Borel transformation with respect to the variable (p + q) 2 in the hadronic and 
the QCD expressions and then equating Eq.(13) with Eq.(18), we obtain the desired sum 
rule for the form factor F s->p (0): 

— -JbF p (0)e t = / <p±(u, u )e T 1 « p'du, (19) 

m b Ju(so) u 

where T is the Borel parameter. 
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It has been shown from the above expression that the distribution wavefunctions ip\\, g\ 
and gj_ disappear in Eq.(19) and the uncertainty due to them is avoided. Thus our sum rule 
for F B ~* P (0) to the twist-3 accuracy is valid, if the corrections from p-meson mass are not 
taken into account. 

Now we will choose some input parameters for the sum rule (19), which are the light-cone 
wavefunction ip±(u, p 2 = 5GeV 2 ), decay constant fs and b-quark pole mass nib- For the same 
reason as in Ref. [6], we can set m& = 4.7, 4.8 and 4.9GeV, accordingly with fs = 165, 120 and 
85MeV. For the light-cone wavefunction (p±(u,p 2 = 5GeV 2 ), taking into account isospin 
symmetry we make use of the same form for p + , p~ , p° and uj: 

<p ± (u, p 2 ) = 6u(l - u)[l + 0.138 x 7.5(£ 2 - 0.2)] (20) 

with f = 2u - 1 in Ref. [10] and 

(p±(u, p 2 ) = 6u(l - «)[1 + 0.077 x 7.5(^ 2 - 0.2) - 0.077 x (39.375^ 4 - 26.25£ 2 + 1.875)] (21) 

with £ = 2u — 1 in Ref. [11] respectively. 

A comparison of the wavefunctions, which are respectively introduced by P.Bali and 
V.M.Braun in Ref. [10] and A.P.Bakulev and S.V.MIkhailov in Ref. [11], is depicted in Fig.l. 
Numerical calculations corresponding to the different wavefunctions will be discussed in the 
following two cases. 

Case a, with the light-cone wavefunction Eq.(20), the best stability in the sum rules for 
F b ^p(0) exists in the window T — 6.0 — 10.0GW 2 , with a threshold s = 32 - 3AGeV 2 . 
Then we obtain the form factors F B ^' 5± (0) = 0.335 ± 0.050, F B ^°^(0) = 0.237 ± 0.035 
and the branching ratio Br(B -> p± 7 ) = (2.71 ± 1.00) x 10~ 6 , Br(B -> (p°,u) + 7) = 
(1.36 ± 0.50) x 10~ 6 . The results are greater than these in Ref. [4] slightly. 

Case b, with the light-cone wavefunction in Ref. [11], the best stability exists in the window 
T = 6.0 — lO.OGeV 2 , with a threshold s = 30 - 33GW 2 . This leads the form factor 
F B ^ ± (0) = 0.272 ± 0.029, F B ^°^(0) = 0.192 ± 0.021 and the branching ratio Br(B -> 
p± 7 ) = (1.79 ±0.61) x 10~ 6 , Br(B -> (p ,cj)+ 7 ) = (0.90 ±0.31) x 10~ 6 . This result accords 
with that F B ^(0) = 0.28 in Ref. [4]. 

5 



The variations of F B ^ P (0) with the threshold Sq and the Borel parameter T are drawn in 
Fig. 2 and Fig.3. Here we would like to emphasize that the form factor F B ^ P (0) is dependent 

2 

on the decay constant /g and the factor m|Exp(^j 1 ) in Eq.(19). This is the reason why 
the large variation of fs with cannot make the form factor F B ~* P (0) a great difference. 
The thresholds sq in our case (a and b) are both lower than that in Ref. [4] because the mass 
of the lowest scalar B meson is smaller than that of the first excited pseudoscalar B meson. 
Compared with the experimental observation in CLEOU, Br(B° — > p 0/ -f) < 3.9 x 10~ 5 , 
Br(B° -> uyy) < 1.3 x 10~ 5 and Br(B~ -> p~i) < 1.1 x 10~ 5 in Ref.[l], our results with the 
two kinds of light-cone wave functions are below the experimental upper limit. 

In summary, we have reanalyzed the rare process B —>■ (p,u) + 7 in the light-cone QCD 
sum rule framework by introducing a chiral current correlator, which is different from that 
in Ref. [4]. In this approach, the chiral-even operator matrix elements, which have not been 
parametrized precisely by light-cone wave functions, disappear for the form factor F B ^ P (0). 
Further, neglecting the corrections from p-meson mass, we have derived the sum rules for 
form factors to twist-3 accuracy. Because of the sensitivity of the numerical results to 
the decay constant /b, we have given precision values of the decay constant fs to leading 
order in QCD, with the b-quark pole mass = 4.7 — A.9GeV. Numerically, the form 
factors, F B ^ P± (0) = 0.335 ± 0.050, F B ^ P °/ W {$) = 0.237 ± 0.035 and the branching ratios 
Br[B -> p ± 7 ) = (2.71 ± 1.00) x lO" 6 , Br(B -> (p°,u) + 7) = (1.36 ± 0.50) x 10~ 6 with 
wavefunctions of P.Ball and V.M.Braun; and the form factors F B ^ P± (0) = 0.272 ± 0.029, 
F B -" P °^(0) = 0.192 ± 0.021 and the branching ratios Br(B -> ^7) = (1.79 ± 0.61) x 10" 6 , 
Br(B -> (p°,w) + 7) = (0.90 ± 0.31) x 10~ 6 with those of A.P.Bakulev and S.V.Mikhailov. 
Both of them are below the experimental upper limit. The different predictions will be tested 
by future precise experiments. 
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Figure captions 

Figure 1: Curves for the different wavefunctions of the p meson. The solid line represents 
the one of P.Ball and V.M.Braun, and the dashed one represents the one of A.P.Bakulev 
and S.V.Mikhailov. 

Figure 2: The form factors F B ~* P (0) corresponding to the threshold s = 33GeV 2 using 
the wavefunction of P.Ball and V.M.Braun. With nib = 4.7,4.8 and A.9GeV respectively, 
the best fits are shown in the window T = 6 — lOGeV 2 . 

Figure 3: The form factors F B ^ P (0) corresponding to the threshold s = 31.5GeV 2 us- 
ing the wavefunction of A.P.Bakulev and S.V.Mikhailov. With m b = 4.7,4.8 and A.9GeV 
respectively, the best fits are shown in the window T = 6 — lOGeV 2 . 
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Fig.l 



